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We develop a three body classical trajectory Monte Carlo (CTMC) method to dicsuss 
the effect of the breakup process on heavy-ion fusion reactions induced by weakly bound 
{Nl '• nuclei. This method follows the classical trajectories of breakup fragments after the 
breakup takes place, and thus provides an unambiguous separation between complete 
^ and incomplete fusion cross sections. Applying this method to the fusion reaction 6 Li + 
209 Bi, we find that there is a significant contribution to the total complete fusion cross 
sections from the process where all the breakup fragments are captured by the target 
nucleus (i.e., the breakup followed by complete fusion). 

1. Introduction 

a 

The fusion reactions of neutron-rich nuclei around the Coulomb barrier provide a good 
opportunity to study the interplay between quantum tunneling and the breakup process. 
An important question here is: does the breakup process influence fusion reactions in a 
similar way as the inelastic process, which leads to subbarrier enhancement [[T] of fusion 
cross sections over predictions for a single barrier? In addressing this question, measure- 
ments with weakly bound stable nuclei, such as 9 Be, 6 Li, and 7 Li, have been proved to 
be useful [0 El II]- These beams are currently much more intense than radioactive beams, 
allowing more precise and extensive experimental studies. Also, these nuclei predomi- 
nantly break up into charge fragments ( 9 Be — > la + n or a+ 5 He, 7 Li —> a + 1, and 6 Li 
— > a + d) , which are more easily detected. This allows a clean separation of the products 
of complete fusion from those of incomplete fusion, where only a part of the projectile is 
captured. 

To describe the effect of breakup on fusion reactions theoretically, a model has to take 
into account the following three different processes: (i) the projectile as a whole is captured 
by the target, (ii) only one of the breakup fragments is captured, (iii) all the breakup 
fragments are absorbed after the breakup takes place near the target nucleus. Process (ii) 
is referred to as an incomplete fusion reaction, while both of processes (i) and (iii) lead 
to complete fusion. Calculations based on the continuum-discretized coupled- channels 
method may account for the total fusion cross section (the sum of all the processes) and/or 
the separate contribution of process (i)[ElE], but it is not easy to compute cross sections 
for process (iii), that is, the breakup followed by complete fusion. In order to model 
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this process, we need to follow the trajectories of the breakup fragments to determine 
whether one or both fragments are captured by the target nucleus. In this contribution, 
we present such a model [|3j|, where many classical trajectories are sampled by the Monte 
Carlo method to compute the total complete fusion cross section. We note that a very 
similar model, called the three-body classical trajectory Monte Carlo (CTMC) method, 
has been developed in the field of atomic and molecular physics, and has been successfully 
applied to ion-atom ionization and charge transfer collisions [|7||B1IH]- 

2. Three-body classical trajectory Monte Carlo (CTMC) model 

Let us assume the following Hamiltonian for a reaction of a projectile nucleus which 
consists of two cluster fragments Pi and P 2 , 

H = + ^ + It + Vn(ri - r2) + VlT(ri - rr) + VMr > - TTh (1) 

where p 1; p 2 , and p T are the momenta of fragment Pi, P2, and of the target, respectively, 
ri, T2-, and vt are the corresponding coordinates. The idea of the three-body classical 
trajectory model is to solve this Hamiltonian classically in the two-dimensional (x, y) plane 
[Ej- We assume that at t=0 the target is at rest at the origin, while the center of mass of 
the projectile is at (x, y) = (P max , b) moving towards the — x direction. Here, P max is some 
large number, while b is the impact parameter. The relative motion between the projectile 
fragments is confined at t — inside the potential Vi 2 at the energy corresponding to the 
breakup Q- value, Qhu, with a random orientation 9 with respect to the x axis. The initial 
conditions thus read, 

xi(0) = P max — — r cos6», 2/i(0) = 6 H — — r sin6>, (2) 

mi + 7712 mi + vri2 

x 2 (0) = P max H ■ r cos6>, 2/2(0) = b ■ r sin^, (3) 

mi + 7B2 mi + m 2 

v x i(0) = -vi ah - ^—p Tcl cos9, v y i(0) = 2~~^ei sin^, (4) 

Vx2{0) = -Vlab + 7T— PrclCOS^, V y 2 (0) = - ^—Pv C \ sin 9, (5) 

ZTYI2 

XT (0) = y T (0)= VxT (0)=v yT (0) = 0, (6) 
where 



^iab = J ■ [Slab - V 1T (r 1T (0)) - V 2T (r 2T (0))}, (7) 

V mi + m 2 



is the initial velocity for the center of mass of the projectile, andp re i = y2/i[<5bu — ^12 (712(0))] 
with /i = mim 2 /(mi + m 2 ). The value of r is restricted to < r < r t , where r t satisfies 
Vi 2 (r t ) = Qbu- 

With these initial conditions, we solve the Newtonian equations to follow the time 
evolution of the trajectories for each particle, Xi(t) , yi(t) , v x i(t) , and v y i(t),(i = 1,2, T). 
At each time, we monitor the relative distance between the particles. We assume that 
the projectile fragment Pj (i = 1,2) is absorbed by the target nucleus when the relative 
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distance between the fragment and the target is smaller than r%r < 1.1(^4/ + A/ ) fm, 
and that the breakup occurs when the relative distance between Pi and Pi is larger than 
the barrier radius for the potential Vyi. By sampling many trajectories for random values 
of r and 9 with the Monte Carlo method, we calculate partial cross sections for a given 
impact parameter b. Total cross sections are then obtained by integrating them over the 
value of b. 
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Figure 1. Examples of classical trajecto- 
ries for the 6 Li+ 209 Bi reaction. Fig. 1(a), 
(b), and (c) are for the impact parame- 
ter 6=2.5, 4.65, and 5.5 fm at E cm =40 
MeV, respectively. The dashed and the 
solid lines correspond to the trajectories 
for the a particle and deuteron, respec- 
tively, while the dotted line is for the cen- 
ter of mass of the projectile fragments. 
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Figure 2. The partial fusion cross sec- 
tions da/db for the 6 Li+ 209 Bi reaction at 
E cm =4:0 MeV obtained with the three- 
body classical trajectory Monte Carlo 
method. 
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Figure 3. The complete fusion (CF) and 
the incomplete fusion (ICF) cross sections 
for the 6 Li+ 209 Bi reaction. The experi- 
mental data are taken from Ref. [Ej- 
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3. Application to 6 Li + 209 Bi reaction 

We now apply the three-body classical trajectory Monte Carlo method to the 6 Li + 
209 Bi system and discuss the effect of breakup of 6 Li into d + a on the fusion reaction 
near the Coulomb barrier. We assume the Woods-Saxon form in calculating V12, Vyn and 
V 2 t, using the same values for the parameters as in Refs. [HUJUI]. We take -R max = 200 
fm. Fig. 1(a) - 1(c) shows sample trajectories at E cm = 40 MeV. These are for b=2.5, 
4.65, and 5.5 fm, respectively. For relatively small values of the impact parameter (Fig. 
1(a)), the projectile nucleus 6 Li is captured by the target before breaking up into a and 
d, leading to complete fusion. For 6=5.5 fm (Fig. 1(c)), only the deuteron is captured 
while the alpha particle escapes from the target. This is an example of incomplete fusion. 
For intermediate values of b, breakup takes place near the target, but both of the breakup 
fragments get absorbed (Fig. 1(b)). This is breakup followed by complete fusion. Figure 
2 shows the partial fusion cross sections for these processes as a function of the impact 
parameter, for E cm =AQ MeV. 

Figure 3 shows fusion cross sections as a function of energy in the center of mass frame. 
The dotted line denotes the sum of complete and incomplete fusion cross sections. The 
dashed line indicates the absorption of the whole projectile from the bound states, while 
the solid line shows the total complete fusion cross sections. The difference between the 
solid and the dashed line thus represents the contribution from breakup followed by fusion 
of both fragments. This makes a significant contribution to complete fusion at energies 
above the barrier. 

4. Summary 

We present a three body classical trajectory Monte Carlo method for fusion of weakly 
bound nuclei. Applying this method to the 6 Li + 209 Bi reaction, we have demonstrated 
that there is a significant contribution of breakup followed by complete fusion to the total 
complete fusion cross sections. This process has not been considered in any coupled- 
channels calculations so far, and the present result suggests that a consistent definition 
of complete fusion is necessary when one compares experimental data with theoretical 
calculations. 
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